home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / EDUCNOMY / ASTROSET.LZH / HELIAC.BAS < prev    next >
BASIC Source File  |  1985-12-20  |  2KB  |  80 lines
  1.  
  2. 10 REM   HELIACAL RISE AND SET
  3. 11 P1=3.14159265: R=P1/180
  4. 12 P2=P1/2: P3=1.5*P1: P4=2*P1
  5. 13 A8=LOG(10): R1=0.017202
  6. 14 INPUT "EXTINCTION";K
  7. 15 INPUT "LIMITING MAG";M0
  8. 16 INPUT "STAR MAG.";M
  9. 17 INPUT "STAR R.A. (HOURS)";A1
  10. 18 INPUT "     DEC. (DEG.) ";D1
  11. 19 INPUT "OBS LATITUDE     ";L
  12. 20 A1=A1*15*R: D1=D1*R: L=L*R
  13. 21 S1=SIN(D1): C1=COS(D1)
  14. 22 S7=SIN(L):  C7=COS(L)
  15. 23 INPUT "RISE OR SET";R$
  16. 24 R$=LEFT$(R$,1)
  17. 25 IF R$="R" THEN F=1
  18. 26 IF R$="S" THEN F=-1
  19. 27 V1=-S7*S1/(C7*C1)
  20. 28 IF V1=0 THEN H1=P1/2
  21. 29 IF V1=0 THEN 32
  22. 30 IF V1<-1 THEN V1=-1
  23. 31 H1=ATN(SQR(1-V1*V1)/V1)
  24. 32 IF V1<0 THEN H1=H1+P1
  25. 33 H1=A1-F*H1
  26. 34 X=-0.2*(M0-7.93+K)
  27. 35 B0=79.4*(10^X-1)^2-K*589
  28. 36 D=80+A1*58.13-60*F
  29. 37 D=D+F: G=D*R1-0.062
  30. 38 PRINT "  TRIAL DAY: ";INT(D)
  31. 39 L0=D*R1-1.406+0.033*SIN(G)
  32. 40 A0=ATN(0.917*SIN(L0)/COS(L0))
  33. 41 D0=0.398*SIN(L0)
  34. 42 D0=ATN(D0/SQR(1-D0*D0))
  35. 43 IF L0<P2 OR L0>P3 THEN 45
  36. 44 A0=A0+P1
  37. 45 IF L0>P3 THEN A0=A0+P4
  38. 46 C0=-S7*SIN(D0)/(C7*COS(D0))
  39. 47 H0=ATN(SQR(1-C0*C0)/C0)
  40. 48 IF C0<0 THEN H0=H0+P1
  41. 49 H0=A0-F*H0
  42. 50 IF F*(H1-H0)>0 THEN 37
  43. 51 N=-100: T=H1-F/30
  44. 52 T=T+F/30
  45. 53 IF F*(T-H0)>0 THEN 37
  46. 54 N0=N: REM  LIMITING MAG.
  47. 55 A9=S7*S1+C7*C1*COS(A1-T)
  48. 56 IF A9<-1 THEN A9=-1
  49. 57 IF A9=0 THEN Z1=P1/2
  50. 58 IF A9=0 THEN 60
  51. 59 Z1=ATN(SQR(1-A9*A9)/A9)
  52. 60 A9=C7*COS(D0)*COS(T-A0)
  53. 61 A9=A9+S7*SIN(D0)
  54. 62 H5=ATN(A9/SQR(1-A9*A9))
  55. 63 A9=C1*COS(D0)*COS(A0-A1)
  56. 64 C2=A9+S1*SIN(D0)
  57. 65 A9=C2-COS(Z1)*SIN(H5)
  58. 66 A9=A9/(SIN(Z1)*COS(H5))
  59. 67 IF A9<-1 THEN A9=-1
  60. 68 T1=ATN(SQR(1-A9*A9)/A9)
  61. 69 A5=0.025*EXP(-11*COS(Z1))
  62. 70 A5=1/(COS(Z1)+A5)
  63. 71 L5=H5*(8.2*Z1+12)+2.86*Z1
  64. 72 L5=4.75-T1*Z1/3+L5
  65. 73 B=B0+(K/0.20)*(10^L5)
  66. 74 IF L5<2.07 THEN B=B0+K*589
  67. 75 C5=4.466E-9: K5=1.258E-6
  68. 76 IF B<1649 THEN C5=1.58E-10
  69. 77 IF B<1649 THEN K5=0.0126
  70. 78 N=C5*(1+SQR(K5*B))^2
  71. 79 N=-16.57-K*A5-2.5*LOG(N)/A8
  72. 80 IF N>N0 THEN 52
  73. 81 IF N0<M THEN 37
  74. 82 IF D>=366 THEN D=D-365.25
  75. 83 PRINT "DAY OF YEAR: ";INT(D)
  76. 84 PRINT "STAR: ";INT(90.5-Z1/R)
  77. 85 PRINT "SUN:  ";INT(H5/R+0.5)
  78. 86 END
  79.  
  80.